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We present an overview of recent results from the CP-PACS computer on the quenched light hadron spectrum 
and an on-going two-flavour full QCD study. We find that our quenched hadron mass results are compatible 
with the mass formulae predicted by quenched chiral perturbation theory, which we adopt in our final analysis. 
Quenched hadron masses in the continuum limit show unambiguous and systematic deviations from experiment. 
For our two-flavour full QCD simulation we present preliminary results on the light hadron spectrum, quark 
masses and the static potential. The question of dynamical sea quark effects in these quantities is discussed. 



1. Introduction 

Calculating the light hadron spectrum from 
first principles has been one of the main objectives 
of lattice QCD since its formulation 0]]. Pioneer- 
ing attempts towards this goal in the quenched 
approximation 0j were followed by an increase of 
understanding and control of various systematic 
uncertainties due to finite lattice spacing, finite 
size and the problem of chiral extrapolations. An 
outstanding effort carried out by the GF11 col- 
laboration|3J led to a quenched spectrum agree- 
ing within 5-10% with experiment. With its his- 
tory quenched QCD can therefore be considered 
as a matured field where an additional increase of 
calculational accuracy will finally tell us how the 
quenched spectrum deviates from experiment. 

Last year the CP-PACS collaboration pre- 
sented first results of an effort in this direction ji]]. 
The simulation has been completed in the mean- 
time and the final results have been presented at 
this conference]^]. 

Progress in the full QCD spectrum study has 
been comparatively slower. Even though full 
QCD simulations have a history almost as long as 
quenched ones[||, the field is still at an exploring 
stage. In particular how in detail the spectrum 
is influenced by the presence of light sea quarks 
has not yet been fully answered [Q. Recently the 
CP-PACS collaboration undertook a further step 
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in the direction of a realistic full QCD simulation. 
At this conference first preliminary results have 
been presented ff-[lj}]- 

In this article we present an overview of the 
two calculations of the CP-PACS collaboration 
described above. In the first part we discuss in de- 
tail our final results for the quenched hadron spec- 
trum, devoting some length to the issue of chiral 
extrapolations. In the second part we switch to 
full QCD where we present our main results, fo- 
cussing on whether effects from sea quarks can be 
seen. 

2. Simulation of Quenched QCD with the 
Standard Wilson Action 

We list the parameters of our quenched simu- 
lation in Table [j], where the partition size of the 
CP-PACS computer employed and the execution 
time for one cycle of configuration generation and 
hadron propagator calculations are also given for 
each lattice size. 

The runs are performed on four lattices cov- 
ering the lattice spacing in the range a -1 = 2- 
4 GeV for extrapolation to the continuum limit. 
The physical lattice size in the spatial direction 
is chosen to be around La rj 3 fm where finite 
size effects are expected to be negligible. Gauge 
configurations are generated with over-relaxation 
and heat bath algorithms mixed in a 4:1 ratio. 

For each value of p quark propagators are cal- 
culated for five quark masses corresponding to 
mp S /m v of about 0.75, 0.7, 0.6, 0.5 and 0.4. The 
first two values are taken to be around the strange 
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Table 1 

Parameters of quenched simulation. Values of 
a -1 are determined from m p obtained by a 
quenched chiral perturbation theory fit to vector 
meson masses. 







[GeV] 



La 
[fm] 



# # 
conf PU 



hrs/ 
conf 



5.90 32 3 x56 1.934(16) 3.26(3) 800 256 3.0 

6.10 40 3 x70 2.540(22) 3.10(3) 600 512 4.8 

6.25 48 3 x84 3.071(34) 3.08(3) 420 1024 6.8 

6.47 64 3 xll2 3.961(79) 3.18(6) 150 2048 15.5 



quark for interpolation, and the last value is our 
attempt to control chiral extrapolations better 
than in previous spectrum studies with the Wil- 
son quark action. We calculate hadron masses 
for equal as well as unequal quark masses so that 
we can determine the complete light hadron spec- 
trum with degenerate u and d quarks and a heav- 
ier strange quark. 

Since Lattice 97 we have increased the number 
of configurations from 91 to 150 at /3 = 6.47. For 
technical details of the simulation we refer to 0] . 

3. Chiral Extrapolation in Quenched QCD 

Extrapolation to the chiral limit is a major 
nontrivial component of light hadron spectrum 
calculations fl2]| . In addition to the difficulty 
that simulations are done far from the chiral 
limit there is the issue of the choice of the func- 
tional form of hadron masses as a function of the 
quark masses. In considering this issue it should 
be asked whether the actual data supports the 
validity of quenched chiral perturbation theory 
(q%PT) (l^ Q , which predicts characteristic sin- 
gularities in hadron masses in the chiral limit. 
One also has to examine how results differ if dif- 
ferent possible fitting functions are employed. In 
the following we examine these questions in de- 
tail. 

3.1. pseudoscalar meson sector 

The dependence of the pseudoscalar (PS) me- 
son mass mps on its two valence quark masses m 
and m s is predicted by q%PT to have the form ]14| ] 

mp s = A(m s + m) {l - S [ln(2mA/A^) 
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Figure 1. Ratio of m PS to quark mass from axial 
Ward identity. Filled symbols are for degenerate 
quark masses, open symbols for non-degenerate 
cases. 
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To test whether our data is consistent with the 
presence of a logarithmic term we plot in Fig. 1 
the ratio of the PS meson mass squared to the 
quark mass. We employ quark masses determined 
from an extended axial current Ward identity 
(AWI), V f A e * t = 2mf Nl aP, avoiding ambiguities 
due to the determination of the critical hopping 
parameter. We see a clear increase of the ratio 
for small quark masses in agreement with ([!]). 

A more quantitative test can be made by form- 
ing the two variables 



= 2- 
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where 7r and r\ are degenerate pseudoscalar 
mesons with an equal quark mass of m or m s , 
and K is the non-degenerate meson. The two 
quantities are related by 



y = 1 + 5 ■ x, 



(4) 



if the 0((m s — m) 2 ) term in (|l|) is ignored. In 
Fig. |^ we plot our mass data in terms of x and 
y for two choices of m s (mps/my — 0.75 or 0.7) 
and three for m (mps/my = 0.6, 0.5 or 0.4). All 
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Figure 2. Pseudoscalar meson test for quenched 
chiral logarithms (see Eqs. (||) and @). 

our data lie within a narrow wedge bounded by 
two lines corresponding to (5 = 0.08-0.12. 

The mass ratio m 2 K /m\ originally proposedjl4| 
has corrections due to both quadratic terms in 
(0) . For this ratio we find that our data does not 
lie on a common line, indicating that the inclusion 
of the term 0((m + m s ) 2 ) is important for the 
description of our data. 

A different test uses a ratio of decay con- 
stants jL{|; the quantity x from (||) and y = 
ffc/ifnlri) are related by y — l—S ■ x/2. As shown 
in Fig. H the data lies again in a wedge, bounded 
by two lines corresponding to (5 = 0.08-0.16. 

Having seen evidence for quenched chiral loga- 
rithms in the above tests, we fit the pseudoscalar 
meson mass with (|l|) , taking into account the cor- 
relation between meson masses at different quark 
masses. We drop the quadratic term proportional 
to C since the mass test @ and (||) suggests its ef- 
fect to be small. The fits are done in terms of l/K 
since it has no statistical error, which means that 
we use the quark mass derived from the vector 
Ward identity (VWI), m^ wl a=(l/K - l/K c )/2. 

In Fig. U we show the degenerate pseudoscalar 
meson mass at = 5.9 as a function of l/K to- 
gether with fit curves from (|l|). Results for a 
quadratic fit in m^ wl a (dashed lines), and the 
AWI quark mass data with a linear fit (dotted 
line) are also plotted. We observe that the AWI 
quark mass is well described by a linear func- 
tion in l/K. This means that a proportional- 
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tions of quark mass, which in turn justifies our 
use of different quark mass definitions for mass 
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Figure 3. Decay constant ratio test. 

tests and fitting. Another important feature is 
a good agreement of the critical hopping param- 
eter K c from the qxPT fit of pseudoscalar me- 
son masses and from a linear fit of AWI quark 
masses, whereas they differ by as much as 17 a if 
a quadratic fit is employed for the pseudoscalar 
meson masses. 

Results for the parameter S obtained from the 
fit are shown in Fig. |5| as function of the lattice 
spacing for degenerate and non-degenerate data 
separately. Values obtained with a linear fit of the 
mass test in Fig. || are also shown (open squares). 
At each value of a, the results of independent fits 
to degenerate and non-degenerate masses agree 
with each other within error bars. The values 
from the mass test are somewhat higher in com- 
parison. Overall however, all the results for 5 are 
mutually in reasonable agreement. 

To summarize our examinations in the pseu- 
doscalar sector, we find the two different direct 
tests to provide a positive evidence for the pres- 
ence of qxPT logarithms. The value of 5 is in the 
range S » 0.10(2). 

3.2. vector meson sector 



For vector meson masses qxPT gives) 16 
m p = my + Ci/ 2 m-n + 2Cim^, 



(5) 



/2 



-(m 7r +m„) + 2 , 
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where C\/2 = —4:irg%5 with g 2 a phenomenological 
coupling constant of the vector meson quenched 
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Figure 4. Chiral extrapolation of the PS meson 
and the AWI quark mass at /3 = 5.9. 
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Figure 5. Parameter S from fit to pseudoscalar 
meson masses (filled symbols) and mass test 
(open symbols). For clarity the data is slightly 
shifted in the x-direction. 



chiral Lagrangian. 

An examination of mass ratios such as (m p — 
m$)l(uiTs — m v ) — C\j2 +0(mps) calculated for 
our data indicates a small negative value in the 
chiral limit. The evidence, however, is less com- 
pelling than for pseudoscalar mesons because a 
chiral extrapolation is needed in order to extract 

Cl/2- 

As an alternative analysis we attempt to fit 
data with the predicted forms (||J|). To examine 
whether our data exhibits the full pattern of de- 
pendence on pseudoscalar masses we employ a si- 
multaneous fit to degenerate and non-degenerate 
masses. Such a combined fit involves a large num- 
ber of data points, for which we find only uncor- 
rected fits to be practically feasible. 

This fit reproduces our vector meson mass data 
well as illustrated in Fig. |[ Clearly visible in 
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Figure 6. Chiral extrapolation of vector meson 
masses at /3 = 5.9. Open symbols in the inset are 
extrapolated values at physical quark mass. 

the fit curve for m p (solid line) is the effect of 
the 0(mps) term, giving rise to an increase close 
to the chiral limit. The magnitude of the effect, 
however, is small, and this holds for all (3. Cor- 
respondingly, small negative values for C-yi 2 are 
found from the fits (see Table |^), consistent with 
indications from mass ratio tests. The average of 
C1/2 over j3 gives Ci/ 2 = —0.071(8), which is ten 
times smaller compared to the phenomenological 
estimate C 1/2 = -^ng^S w -0.71 (for g 2 = 0.75 
and 6 = 0.1). 

3.3. baryon sector 

For the baryons we find a situation for mass 
ratio tests less conclusive than for vector mesons. 
Therefore we follow the same strategy and fit 
masses of different degeneracies simultaneously. 
We assume qx?T mass formulae i.e. for 

octet baryons 



(AF 2 w uu - A(D - F)Fw u 
+ (D — F) 2 w ss ) /2 
4b F m 2 uu + 2(b D - b F )m 2 ss 
(2D 2 /3~2F 2 )v uu 
(2D 2 /3-4DF + 2F 2 )v us 
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Figure 7. Chiral extrapolation of baryon masses 
at (3 — 5.9. Open symbols in insets are extrapo- 
lated values at the physical quark mass. 

- 2(4D/3 - 2F)(D/Z + F)w us ) /2 
+ 4(26^/3 - b F )m 2 uu - 2(b D /3 + b F )m 2 ss 
+ {2D 2 /9-8DF/3 + 2F 2 )v uu 
+ {10D 2 /9-4DF/3-2F 2 )v us , 

and for decuplet baryons 
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where m uu and m ss are n and r\ meson masses 
and w us = -2TTS(m 3 ss -ml u )/(m 2 s -ml u ), v us = 
(m^ s /167r/ 2 ) and / is the pion decay constant. 

These forms are slightly simplified compared 
to the full formula in order to get stable fits. For 
decuplet baryons all leading linear terms in mps 
and all quadratic terms are included, but not the 
cubic terms 0(mp s ). For octet baryons, on the 
other hand, a bending in the nucleon mass data 
necessitates a retention of such terms; we keep 
those arising from the octet-octet couplings F and 



Table 2 

Coefficient C\/2 of the mps term for vector me- 
son and baryons. 



P 

N 

A 



= 5.90 = 6.10 /3 = 6.25 j3 = 6.47 



-0.058(27) -0.075(28) -0.065(35) -0.155(72) 
-0.126(15) -0.120(15) -0.116(13) -0.087(30) 
-0.169(33) -0.145(39) -0.101(40) -0.184(104) 



D, but not those due to the decuplet-octet cou- 
pling C. We also note that / = 93 MeV and 5 = 0.1 
are fixed in the fits. 

The result of the mass fits to octet and de- 
cuplet baryons is shown at f3 = 5.9 in Fig. [t|. 
The increase of masses close to the chiral limit 
as a quenched artifact is visible in the fit curves 
for the nucleon, E, A and £*, but is small com- 
pared to the error after extrapolation. For the nu- 
cleon the coefficient of the 0(mps) term C1/2 = 
— 3n(D — 3F) 2 5/2 obtained from the fits has an 
average value -0.118(4) (see also Table |). This 
is considerably smaller than a phenomenological 
estimate —0.27 for which 6 = 0.1, and F = 0.5, 
Z? = 0.75 estimated from experiment fl8f| are used. 
The contribution of this term to the nucleon mass 
is therefore about 2% or 15 MeV. A similar value 
is found for A. 

Summarizing, our data does not provide di- 
rect evidence for the validity of qxPT in vec- 
tor meson and baryon sectors. However, qxPT 
formulae describe the overall pattern of the data 
well, with the fit parameters taking reasonably 
constant values at different values of f3. Our 
data is therefore compatible with the presence 
of 0(mps) terms. Their magnitudes are smaller 
than expected, amounting to about 2% of vector 
and baryon masses. 

3.4. polynomial chiral fits 

In order to investigate the effect of choosing 
different fitting functions we have also made chi- 
ral extrapolations employing polynomials in 1 / K 
(cubic for N, quadratic otherwise). Hadrons 
of different degeneracies are fitted independently 
and correlations between different quark masses 
are taken into account. The results of these fits 
are already plotted for each channel in Figs. |J, 
H and by dashed lines. The continuum ex- 
trapolations with qxPT and polynomial fits are 
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Figure 8. Continuum extrapolation of quenched 
hadron masses with Mk as input for qxPT and 
polynomial chiral fits. 

compared in Fig. || for meson masses and typical 
baryon masses. 

The two types of fits do lead to different results 
after chiral extrapolation at each f3, which can be 
as significant as 5er in the largest case. In absolute 
magnitude, however, this difference corresponds 
to 2-3% at most. Furthermore, after continuum 
extrapolation, the results from polynomial fits do 
not differ beyond 1.5% from the values obtained 
with qxPT fits, which is a 1.5er effect. 

We conclude that qxPT and polynomial chi- 
ral extrapolations, while possibly leading to siz- 
ably different continuum values in principle, ac- 
tually yield results which do not differ beyond a 
2% level. 

4. Quenched Light Hadron Spectrum 

Based on the results of analyses described in 
Sec. 3 we choose chiral extrapolations using qxPT 
formulae for the calculation of our final spec- 
trum results. The physical point for the degen- 
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Figure 9. Final results of the quenched light 
hadron spectrum in the continuum limit. Ex- 
perimental values and results from GF11 J3| are 
plotted for comparison. 



erate u and d quarks is fixed by M n (135.0) and 
M p (768.4), and we use M^(497.7) or M (1O19.4) 
for the s quark. The lattice scale is set with M p . 

We observe in Fig. || that the continuum ex- 
trapolation is well described by a linear function 
in the lattice spacing. Indeed, if we parameterize 
the a-dependence as m = m (l + aa), we find 
the typical scale a to be a ~ 0.2 GeV. Hence 
higher order terms are expected to be small, e.g., 
(aa) 2 ~0.01 at a = 0.5GeV _1 . We thus employ a 
linear continuum extrapolation. We present our 
final result for the quenched light hadron spec- 
trum in Fig. ^. 

This figure shows that the quenched light 
hadron spectrum systematically deviates from ex- 
periment. If one uses Mk as input to fix the s 
quark mass, the mass splitting between K and 
K* (</>) in the meson sector is underestimated by 
9.5% (12%), which is a 4.3ct (4.9a) effect. In the 
baryon sector, the octet masses are smaller, by 
7.6% (4.5er) on the average, and the decuplet mass 
splittings are 30% (I.Ict for individual level spac- 
ings, but 3.2u for the A — 57 splitting) short of 
the experimental splitting. With as input, 
the K* meson mass is in agreement with experi- 
ment, and the discrepancies in the baryon masses 
are reduced. However, the K — K* meson hyper- 
fine splitting is still smaller by 16% which has a 
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Figure 10. Continuum extrapolation of light 
quark masses in the MS scheme at /i = 2 GeV. 

statistical significance of 6.1a. 

These conclusions remain unchanged compared 
to those presented at Lattice 97 [f|, where we em- 
ployed a linear chiral extrapolation in 1 / K (cubic 
for TV and quadratic for A), except in numerical 
details for individual channels. The most conspic- 
uous change is the decrease of the nucleon mass 
away from the experimental value, and that of the 
A mass toward it. Strange baryon masses with 
Mk as input have also decreased. The changes, 
however, are at most 1.5a of the error of our final 
results, and do not affect the overall pattern of 
deviation between the quenched and experimen- 
tal spectrum. 

5. Quark Masses in Quenched QCD 

The use of qxPT fits to pseudoscalar me- 
son masses has important consequences for VWI 
quark masses. As already noted with Fig. || the 
critical hopping parameter is different for polyno- 
mial and qxPT fits of pseudoscalar meson masses. 
This translates into a different estimate of m^ wl , 
influencing the average u and d quark mass most. 

In Fig. |l(^ we show the quark mass results as 
a function of the lattice spacing. The conclu- 
sion, reported at Lattice 97, that the VWI and 
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Figure 11. Pseudoscalar meson decay constants. 



AWI quark masses converge in the continuum 
limit remains unchanged also when q%PT fits are 
used. Another point to note is that the continuum 
result of the VWI quark masses obtained with 
quadratic polynomial chiral fits of pseudoscalar 
meson masses are consistent within error bars 
with those of qxPT fits, which is similar to the sit- 
uation with hadron masses discussed in Sec. 3.4. 



For our final values of the quark masses we 
make a combined linear continuum extrapolation 
of the VWI and AWI quark masses from qxPT 
fits together, restricting them to a common value 
in the continuum limit (see Fig. 3 in (5)). We 
obtain m ud = 4.6(2) MeV and m s = 115(2) MeV 
[M K input) or m s = 143(6) MeV (M^ input) in 
the MS scheme at /1 = 2 GeV. The value of m u d 
is 12% larger than the estimate reported at Lat- 
tice 97, where linear chiral extrapolations of pseu- 
doscalar meson masses have been used. 



6. Pseudoscalar Meson Decay Constants 

PS meson decay constants are determined from 
the local axial current. Because the quark mass 
dependence of degenerate decay constants has not 
yet been determined within q%PT, we employ a 
quadratic polynomial chiral extrapolation for fa 
and a linear extrapolation for fx- The results 
as a function of the lattice spacing are shown in 
Fig. !| We obtain fa = 120(6) MeV and f K = 
139(4) MeV in the continuum limit, which are 10 
and 15% smaller than experiment, respectively. 
The ratio fx/ fa — 1 = 0.156(29) is also smaller 
than the experimental value 0.212. 



8 



7. Full QCD Simulation 

Having observed systematic deviations of the 
quenched hadron spectrum from experiment, the 
natural question to ask is whether the deviation 
can be explained as a consequence of the effect of 
dynamical quarks. Due to the large increase of 
computer time required for simulating full QCD 
compared to the quenched approximation, we are 
forced to work at coarse lattice spacings where 
a -1 < 2 GeV, and are therefore led to employ 
improved actions. 

In an initial study Jl9[ ] we compared the scaling 
behaviour of various combinations of improved 
gauge and quark actions in full QCD at lattice 
spacings a -1 w 1 GeV. Based on the outcome of 
this comparison we decided to employ for gauge 
fields the RG improved action consisting of pla- 
quette and 1x2 rectangular Wilson loop proposed 
by Iwasaki [EC|. For quarks we use the SW 
clover action f2lfl with mean-field improved value 
of c sw - c SW = P~ 3/4: with P = 1 - 0.8412 • 
the one-loop value of the plaquette. 

Simulations are carried out with two flavours 
of sea quarks, to be identified with the degener- 
ate u and d quarks, while the strange and heav- 
ier quarks are treated in the quenched approxi- 
mation. In Table || we give an overview of our 
run parameters. Three lattice spacings in the 
range a -1 ss 1-2 GeV are used for continuum ex- 
trapolation. Lattice sizes are chosen so that the 
spatial size is kept around La w 2.4 fm (except 
for (3 — 2.2 where it turned out to be smaller). 
For chiral extrapolations, simulations are per- 
formed at four sea quark masses corresponding 
to mps/mv w 0.8, 0.75, 0.7 and 0.6. 

Configurations are generated with the HMC al- 
gorithm. To speed up our code we have imple- 
mented several improvements. For the integra- 
tion of molecular dynamics equations we employ 
the scheme 

Tri-d^Tui^Tpiil-bWTu^Tp^dT), 

where Tp(dr) moves the gauge field and Tu(dr) 
the conjugate momenta by a step dr. This scheme 
has errors of the same order as the standard 
leapfrog but the main contribution to the error is 
removed by the choice b = (3 — \/3) /3 . For our 



Table 3 



Parameters of full QCD simulation. 



O i 70 n — i rp n vl 
bize a v j 

#PU La [fm] 

P 

csw 


Ksea, Trips #traj time 
/my /traj 
[hrs] 


i/ xz4 o.yi((ioj 
64 2.58(3) 
1.80 
1.60 


.1409 .8060(/J 6250 0.11 
.1430 .753(1) 5000 0.16 
.1445 .696(2) 7000 0.27 
.1464 .548(4) 5250 0.93 


16 a x32 1.288(15) 
256 2.45(3) 
1.95 
1.53 


.1375 .8048(9) 7000 0.11 
.1390 .751(1) 7000 0.16 
.1400 .688(1) 7000 0.26 
.1410 .586(3) 5000 0.85 


24 3 x48 2.45(9) 
512 1.93(7) 
2.20 
1.44 


.1351 .800(2) 1250 0.44 
.1358 .752(3) 1350 0.70 
.1363 .702(3) 1610 0.99 
.1368 .637(6) 1265 1.98 


heaviest quarks dr can be taken about a factor 3 



larger than for leapfrog, leading to a 30% gain in 
computer time. The gain decreases for decreas- 
ing quark masses and is lost at our smallest quark 
mass. For the inversion of the quark matrix we 
use the even/odd preconditioned BiCGStab algo- 
rithm p3fl . Its performance compared to MR in- 
creases with decreasing quark mass. We observe a 
gain of about 40% at m P s/mv~0-7. CPU times 
needed for the creation of one unit of trajectory 
and numbers of processors are listed in Table ||. 

The light hadron spectrum is calculated at each 
sea quark mass with valence quark masses corre- 
sponding to mps/mv ~ 0.8, 0.75, 0.7, 0.6 and 
0.5. Quark propagators are calculated with ex- 
ponentially smeared sources and local sinks. We 
combine quark propagators to form mesons of the 
form SS, SV and VV and baryons of the form 
SSS, SSV, SVV and VVV where S stands for 
a valence quark with K va i = K sea and V for a 
valence quark with K va i ^ K sea . These combi- 
nations include all unequal quark mass combina- 
tions allowed for degenerate u and d sea quarks 
and valence s quark. 

We measure hadron propagators at every 5th 
trajectory. Hadron masses are extracted by un- 
corrected fits. Errors are determined with the 
binned jackknife method. The error reaches an 
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approximate plateau after a bin size of at most 
10 configurations (50 HMC trajectories), which 
we choose for error estimates of all our runs. 



8. Full QCD Light Hadron Spectrum 

Figure [l^ shows meson masses at /3 = 1.95 
as a function of the average hopping parameter 
1/K aver = (1/K val(1) + l/K val{2) )/2 of the two 
valence quarks in the meson. Partially quenched 
data (VV and SV) lie along clearly distinct lines 
for each K sea . Each of these mass results are 
almost linear for both pseudoscalar and vector 
mesons, but their slope shows a variation with 
K sea . With decreasing sea quark mass we ob- 
serve a decreasing slope for pseudoscalar mesons 
and an increasing slope for vector mesons. While 
the pseudoscalar meson mass in the S S channel 
is almost linear in X/K sea , the vector meson mass 
exhibits a clear concave curvature. 

The picture for decuplet baryons looks very 
similar to the one for vector mesons, whereas for 
octet baryons partially quenched data of different 
degeneracies lie on clearly different lines. 

For chiral extrapolations we fit at each value 
of j3 all the data of a given channel together in 
a combined fit with a general quadratic ansatz 
in 1/K sea and 1/K va i. For pseudoscalar mesons 
consisting of valence quarks "1" and "2" , the fit 



ansatz is 

m 2 PS = B s m sea + B v m va i + C s m 2 sea + C v m 2 val 
+C sv m sea m va i + Ci2m va i(X)m va i(2) (10) 

where bare quark masses are defined as 

m sea = {l/K sea -l/K c )/2 (11) 
m val{{) = (l/K val(l) -l/K c )/2 (12) 

and m va i is the average valence quark mass. The 
valence-valence cross term C\i is introduced to 
take into account mass differences of mesons at 
the same average valence quark mass. For vec- 
tor mesons and decuplet baryons such a term is 
dropped because we observe them to be well de- 
scribed as a function of the sea quark mass and 
the average valence quark mass only. Curves of 
the combined fit are illustrated in Fig. [l^ for me- 
son masses. 

For octet baryons we employ a combined 
quadratic fit for E-like and A-like baryons simul- 
taneously. Linear terms are taken in a form in- 
spired by chiral perturbation theory and general 
quadratic terms of individual rn va \u\ are added. 
The resulting ansatz has 12 parameters which fits 
the 84 octet baryon mass data per /3 value well. 

Using the results from the combined quadratic 
fits, the physical light quark hopping parameter 
K u d and the lattice scale a -1 are fixed with the 
experimental value of and M p as input. The 
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Figure 13. Hadron masses in full QCD as func- 
tion of the lattice spacing. Mk is used as input. 
At the smallest lattice spacing A, E, E* and 2* 
masses arc shifted in x-direction for better visibil- 
ity. For comparison quenched QCD results with 
the Wilson quark action are included. 

valence strange quark hopping parameter K s (at 
Ksea = K ud ) is determined from either M K or 
M^. The light hadron spectrum is calculated by 
setting K sea to K ud , and K val to either K ud or K s 
depending on the quark content of the particle. 
The resulting hadron masses as a function of the 
lattice spacing are shown in Fig. [l^. In the next 
section we discuss these results, focussing on what 
extent dynamical sea quark effects can be seen in 
comparison to the results from our simulation of 
quenched QCD with the standard action, plotted 
in Fig. 03 by open symbols. 



9. Sea Quark Effects 

9.1. effects on the light hadron spectrum 

An interesting result which stands out in 
Fig. |l^ is that the K* and <fi masses from the 
full QCD simulation lie higher than those of the 
quenched results at finite lattice spacing, and that 
they extrapolate to values in the continuum limit 
considerably closer to experiment. If one makes 
a linear continuum extrapolation as shown by 
dashed lines, the difference from the experimen- 
tal value is decreased by about 60% compared to 
quenched QCD. The expected leading scaling vio- 
lation for the tadpole-improved SW quark action 
is 0(g 2 (a)a). We also try such an extrapolation, 
which yields similar results. 

The origin of the remaining difference from ex- 
periment is an interesting issue. One possible rea- 
son is that sea quark masses explored in the sim- 
ulation (0.8 > mps/mv ^ 0.6) are still relatively 
heavy so that sea quark effects arc only partly 
operative. We should also remember that our 
simulation treats only two quark flavours as dy- 
namical. Therefore, another possibility is that at 
least part of the difference reflects the quenched 
treatment of the strange quark. 

It is interesting to see in detail how the dif- 
ference between quenched and full QCD arises 
as a function of the sea quark mass. The J pa- 
rameter J = my dmv / dm 2 P s at my / 'rap s = 1.8 
p4| ] is a measure of the meson hyperfine split- 
ting around the strange quark. With the strange 
quark treated in the partially quenched approx- 
imation the derivative is evaluated at fixed sea 
quark mass by varying the valence quark mass 
only. Decreasing the sea quark from a large value 
down to the light quark mass interpolates be- 
tween quenched and two-flavour full QCD. 

In Fig. |l4] we plot J as a function of the sea 
quark mass (in the plot the mass squared of the 
pseudoscalar meson consisting of two sea quarks 
is used). We find the values to be approxi- 
mately constant for moderate or large sea quark 
masses and roughly consistent with the quenched 
value 0.346(22) in the limit of infinitely heavy sea 
quarks. Only when the sea quark mass becomes 
close to the chiral limit the value of J rises to- 
ward the experimental value of 0.48, resulting in 
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Figure 14. J parameter at different sea quark 
masses and lattice spacings. The experimental 
value is denoted in the chiral limit with a star. 

an increase of the hyperfine splitting. 

Sea quark effects on baryons are less clear. 
While scaling violation appears small for baryon 
mass results, the errors are still too large to reli- 
ably discuss continuum extrapolations and how 
the two-flavour results compare with those of 
quenched QCD. Furthermore there are two points 
that have to be kept in mind: (i) The spatial 
lattice size at the smallest lattice spacing turned 
out to be 1.9 fm compared to 2.5-2.6 fm at the 
two coarser lattice spacings. These are smaller 
than the size of 3 fm for our quenched simula- 
tion, and raises the question of finite-size effects, 
(ii) The smallest valence quark mass for full QCD 
corresponds to mps/mv ~ 0.5 rather than rs 0.4 
for quenched QCD, possibly leading to systematic 
errors in chiral extrapolations of baryon masses. 
Further work matching these points is needed to 
clarify sea quark effects in the baryon mass spec- 
trum. 

9.2. static quark potential 

In full QCD the string between a static quark 
and antiquark can break and one would expect 
to observe a flattening of the potential at large 
distances. 

We extract the value of the potential between 
static quarks V(R) from smeared Wilson loops. 
At /3 = 1.95 we also measure the mass of the 
static-light pseudoscalar meson itisl with the 
light quark mass equal to the sea quark mass. 
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Figure 15. Static quark potential and static-light 
meson mass at j3 = 1.95 and the lightest sea quark 
mass. 

In Fig. |l5| these two quantities are compared for 
the lightest sea quark mass at mps/my ~ 0.6. 
The linear increase of the potential continues be- 
yond the crossing point with 2msL at around 1 
fm. A similar situation is found also at other sea 
quark masses. We therefore see no indication for 
string breaking, similar to the findings of previous 
lattice simulations 0. 

A possible explanation of this result is that the 
Wilson loop does not have a sufficient overlap 
with the state of a broken string. In Q we argue, 
that in order to find the true ground state one 
would have to extract the potential from Wilson 
loops at a large time separation. This is diffi- 
cult in practice because of the rapidly vanishing 
signal in the noise. One would have to use (an) 
operator (s) with a better overlap with a static- 
light meson pair state. This is also the essence of 
other recent proposals p5[. 

A quantity in which clear sea quark effects can 
be seen is the Coulomb coefficient a. Here the full 
QCD results are larger than quenched, and they 
increase with decreasing sea quark mass. The 
magnitude of the shift from quenched to full QCD 
is consistent with expectations for Nj — 2 || . 

10. Light Quark Masses in Full QCD 

It has been known for some time that light 
quark masses calculated in two-flavour full QCD 
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take significantly smaller values compared to 
those in quenched QCD|26jj. Most recently 
the SESAM collaborationf27|| reported m u d = 
2.7(2) MeV from their simulation with the stan- 
dard action at a" 1 « 2.5 GeV (j3 = 5.6). They 
furthermore pointed out that the ratio m s /m u d 
takes a large value ~ 50 in contrast to the value 
~ 25 commonly quoted from chiral perturbation 
theory. 

These results have been obtained with the VWI 
definition of the quark mass m^ WI = Z m (l/K q — 
l/K crit )/2, where K^ ea , the hopping parameter 
where the PS meson mass vanishes when valence 
and sea quarks have the same mass, is substituted 
for K crlt . 

It has been suggested Q that the use of the 
partially quenched critical value defined by 

mps(K va i) = with K sea = K u d fixed, resolves 
the problems described above. In two-flavour 
QCD however, the choice K s c ea is the most natu- 
ral one for the average u and d quark mass m u d, 
since the sea quarks are identified with the degen- 
erate u and d quark. Furthermore, K^ ea and Kf?® 
coincide in the continuum limit. If the difference 
between them vanishes sufficiently fast when ap- 
proaching the continuum limit, then also the two 
m^ wl will agree. Indeed, calculating m^ WI in 
perturbation theory, we find that the difference 
between the two quark mass definitions vanishes 
as 0(g i a log a). 

Figure |l6| a shows our results for m u d in the MS 
scheme at fj, — 2 GeV, calculated with tadpole- 
improved one-loop values of renormalization fac- 
tors |ll|, as a function of the lattice spacing. The 
two definitions of m^ WI lead to significantly dif- 
ferent values at finite lattice spacings but as 
expected they converge towards the continuum 
limit. In addition the AWI quark mass also ex- 
trapolates to a similar value. The value reached 
is 2.5-3 MeV, which is about 40% smaller than 
the result obtained in quenched QCD. A theo- 
retical argument for a decrease of quark masses 
by unquenching has been previously made from 
different running of the f3 function |29||. 

Results for the strange quark mass are shown 
in Fig. |l6|b. The difference between various quark 
mass definitions is much smaller than for the 
light quark. The main uncertainty stems from 
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Figure 16. Quark mass in full QCD with various 
definitions. 

the choice of the experimental input to fix K s . 
We find values in the continuum limit of about 
m s « 70 MeV (M K input) or 80 MeV (M in- 
put). The discrepancy between these two values 
is smaller than in quenched QCD. The remain- 
ing difference might be due to the quenching of 
the strange quark or that the sea quarks are too 
heavy. 

In Fig. [l6|c we show the ratio of strange to 
light quark masses as function of the lattice spac- 
ing. While for VWI masses with the choice 
K cr it — K & c ea the value is ~ 45 at the largest lattice 
spacing, it decreases toward the continuum limit 
and all definitions converge to m s /m u d ~ 25, con- 
sistent with xPT results. 
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11. Conclusions 

In our extensive study of the quenched hadron 
mass spectrum, we found that the mass results 
show an overall consistency with qxPT predic- 
tions, which led us to employ qxPT mass formu- 
lae for chiral extrapolations in our final analysis. 
The spectrum in the continuum limit deviate sys- 
tematically and unambiguously from experiment, 
on a level of about 10% (5a). The terms spe- 
cific to qxPT, however, have a small magnitude 
of 2-3% of hadron masses, and are hence respon- 
sible only for a small part of the deviation. While 
chiral extrapolations differ, our main conclusion 
remains unchanged since Lattice 97. 

Our study of full QCD is in a less advanced 
stage, attempting to systematically approach a 
realistic QCD simulation. Still, already several 
intriguing results are found. The discrepancy 
in the meson hyperfine splitting is reduced by 
about 60% and light quark masses are found to 
be around 40% smaller than in quenched QCD. 
In the future it will be interesting to look for sea 
quark effects also in other areas such as topology 
and heavy quark physics. Work in these direc- 
tions has already started. 

I would like to thank the members of the CP- 
PACS collaboration for their continuous help and 
support during the preparation of this overview. 
This work is in part supported by the Grant-in- 
Aid No. 09304029 of the Ministry of Education. 
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